Nonequilibrium stationary states with ratchet effect 
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An ensemble of particles in thermal equilibrium at temperature T, modeled by Nose-Hoover dy- 
namics, moves on a triangular lattice of oriented semi-disk elastic scatterers. Despite the scatterer 
asymmetry a directed transport is clearly ruled out by the second law of thermodynamics. Introduc- 
tion of a polarized zero mean monochromatic field creates a directed stationary flow with nontrivial 
dependence on temperature and field parameters. We give a theoretical estimate of directed current 
induced by a microwave field in an antidot superlattice in semiconductor heterostructures. 
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According to the second law of thermodynamics there 
is no stationary directed transport in spatially periodic 
asymmetric systems in thermal equilibrium |l|, |2| . How- 
ever, a time periodic parameter variation may drive such 
a system out of equilibrium leading to the emergence of 
stationary transport whose direction depends non triv- 
ially on parameters. Such directed transport appears in 
systems with noise, fluctuations and dissipation and is 
now called Brownian motor or ratchet (see e.g. reviews 
01 ) The ratchet effect has a generic nature and 
it has been observed in various physical systems includ- 
ing semiconductor heterostructures cold atoms in a 
laser field 7] , vortices in superconductors |3, IE LLQ| and 
macroporous silicon membranes under pressure oscilla- 
tions . It has also important applications in biological 
systems as discussed in |j, |l2| . 

In spite of a great recent interest to ratchets the 
theoretical research is mainly concentrated on one- 
dimensional models (see e.g. |g). Also, since the ratchet 
behavior is usually rather complex, an overdamped limit 
is used very often to obtain analytical parameter depen- 
dence even if in this regime a directed transport is absent 
for ac zero mean force 5] . To understand in a better way 
the global properties of ratchets and their dependence 
on such important physical parameters as temperature 
T and driving strength f , we analyze here a generic case 
when ac-driving affects a Maxwell thermostat ensemble 
of noninteracting particles moving in an asymmetric two- 
dimensional (2D) periodic structure. This structure is 
composed of triangular 2D-lattice of rigid semi-disks of 
radius as shown in Fig. 1 (insert). The distance R 
between disk centers is fixed to be R = 2r^ and we 
assume that collisions with semi-disks are elastic. Free 
particle motion between semi-disks is affected by a polar- 
ized monochromatic force f = /(cos#, sin#) coscjt with 
frequency lj, strength / and polarization angle 6 to x- 
axis. It is also assumed that particles are in thermal 
equilibrium and at / = their velocities are given by 



the Maxwell distribution at temperature T. Fig. 1 shows 
that in this system ac-force generates stationary directed 
transport. In numerical computations we put r^, particle 
mass m, unit of time and Boltzmann's constant k to be 
equal to unity. 




FIG. 1: (color online) Density distribution averaged over the 
time interval < t < 5 • 10 5 and obtained from dynamics 
of 200 particles given by the Nose-Hoover equations at ther- 
mostat temperature T = 24. The region of distribution is 
x = [-2050,150], y = [-300,1900]. Initially particles are 
placed at x = y = (cross) with random velocities. Den- 
sity is proportional to color changing from zero (red/black) 
to maximum (yellow/ white). The parameters of driving force 
are / = 16, a; = 1.5 and = tt/8. The relaxation time scale 
of the thermostat is r — \/50- Insert shows one trajectory on 
small scale moving between semi-disks. 



In order to put particles in thermal equilibrium we 
choose the elegant method of the Nose-Hoover thermo- 





FIG. 2: (color online) Steady state distribution in 2D momen- 
tum plane (p x ,p y ), density is proportional to color changing 
from zero (blue/black) to maximum (rose- violet/gray) . Left: 
the Maxwell distribution at temperature of Fig. 1; right: dis- 
tribution obtained numerically from the Nose-Hoover thermo- 
stat for the case of Fig. 1. 

stat (see e.g. [ill Q, [H3 and Refs. therein). In this 
method the motion of a particle is affected by an effective 
friction 7 which keeps the average kinetic energy (p 2 /2) 
equal to a given thermostat temperature T. In this way 
the dynamics of particle is described by the equations: 

q = p/m, p = F- 7 p, 7 = [p 2 /(2mT) - l]/r 2 (1) 

where q, p are particle coordinate and momentum, F is 
a sum of ac-force and force of elastic collisions with semi- 
disks, and r is the time scale of relaxation to equilibrium. 

It is known that the Nose-Hoover thermostat works 
well only if the dynamics is sufficiently chaotic 0, Q, 
Hij . In some cases, e.g. for the Galton board, the Nose- 
Hoover thermostat gives noticeable deviations from the 
Maxwell distribution . To check that in our case this 
method really gives a thermal equilibrium we analyze 
the steady state distribution in the momentum space ob- 
tained by numerical Runge-Kutta integration of Eqs.Q. 
Our results show that a small ac-force is needed to make 
chaotic dynamics between semi-disks more homogeneous 
and to produce a stable Maxwell thermal equilibrium 
which is not sensitive to variation of relaxation rate 1/r 
(Fig. 3, insert). At large force the numerical data show 
that 2D steady state in the momentum space is still close 
to the Maxwell distribution (Fig. 2) even if the ac-driving 
produces a clear ratchet effect shown in Fig. 1. The de- 
pendence of steady state on temperature closely follows 
the Maxwell distribution in momentum space p = |p| as 
shown in Fig. 3. Thus we may conclude that the dynam- 
ics given by the Nose- Hoover equations allows efficiently 
investigate the effects of ac-driving on particles in ther- 
mal equilibrium. The numerical data show that this driv- 
ing generates a strong ratchet effect (Fig. 1) with directed 
transport wich depends on temperature and parameters 
of driving force. 

To understand the properties of this directed transport 
we first analyze the dependence of averaged friction (7) 




FIG. 3: (color online) Thermal distribution p in momentum 
p = |p| for different values of temperature T = 8, 16, 40 (from 
narrow to broad distribution, respectively, shown by curves 
and symbols). For each temperature value the curve gives 
the Maxwell distribution and the symbols show the numerical 
data for the Nose-Hoover thermostat (r = a/50) in presence 
of ac-driving with parameters of Fig. 1. Insert shows the 
stability of Nose-Hoover thermostat at small ac-force (/ = 
0.5, 00 = 1.5, = 7r/8) with respect to variation of relaxation 
time t 2 = 1 (circles), 50 (squares), 100 (diamonds) at fixed 
temperature T — 24; curve shows the Maxwell distribution. 
All numerical data are obtained from one long trajectory with 
t<5-10 5 . 

on driving strength / and temperature T. The value 
of (7) is obtained by averaging over long time interval 
during numerical integration of Eqs.Q for one trajectory. 
We also checked that averaging over a few trajectories 
gives the same result. The data are shown in Fig. 4. 
They are well described by a global scaling given by 

( 7 > = Cram-^f/T^ . (2) 

Small deviations seen at low T appear because of strong 
driving force which starts to modify significantly the par- 
ticle velocity distribution in this regime. The numerical 
constant C is only weakly dependent on and uo chang- 
ing by 50% to 30% when changes from to tt/2 and u 
changes by a factor 10, respectively. The dependence Q 
clearly tells that in presence of driving force the thermo- 
stat creates an effective friction force = — (7)p acting 
on particle propagation with an effective friction constant 
(7). Surprisingly, this friction coefficient varies with / 
and T according to Eq.Q but in a large range remains 
independent of the relaxation time r (Fig. 4 insert). We 
note that the particle dynamics in absence of thermostat 
but in presence of friction force f/ = — 7P with constant 
friction coefficient 7 has been analyzed in 16] where it 
was shown that ac-force generates a directed transport 
on semi-disk lattice. 

To understand the origin of the dependence (0 we put 
forward the following heuristic arguments. The driving 
force gives a diffusive energy growth during a dissipative 
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FIG. 4: (color online) Dependence of rescaled average fric- 
tion coefficient (7) / f 2 on temperature T for / =4 (circles), 
8 (squares), 16 (diamonds) 32 (triangles) (top to bottom) at 
fixed uj — 1.5, # = 0, r 2 = 50. The straight line shows de- 
pendence with C = 0.02. Insert shows that (7) is robust 
against variation of r 2 in the interval [1,100] (with step 1); 
data are shown for / = 16, uj = 1.5, = 0, T = 8 (top curve) 
and 24 (bottom curve). 



time scale I/7 so that 

(AE) 2 ~D E ft, D E ~f 2 vl, (3) 

where the diffusion rate in energy is De ~ E 2 r c ~ f 2 vl 
and the mean-free path I ~ R ~ r& ~ 1 determines 
the collision time r c = l/v. In the Maxwell equilib- 
rium the particle velocity is v ~ (Tim) 1 1 2 and the fact 
that the driving force does not modify the velocity dis- 
tribution implies that AE ~ T so that the diffusive 
growth is stopped by effective friction 7 ~ D E /T 2 ~ 
r^m _1 / 2 / 2 /T 3 / 2 in agreement with (|2J). In fact there is 
a close relation to results ^(| where the thermostat is 
absent but a friction force ff = — 7P with constant 7 
affects particle dynamics. In that case ac-driving force 
heats a particle up to energy E ~ (raf 2 /m 1 / 2 ^) 2 / 3 while 
in presence of thermostat the energy is fixed by tempera- 
ture T ~ E that imposes a convergence to the stationary 
state with effective friction given by (0) 

The dependence of average velocity Vf of the ratchet 
flow on (7) for various values of driving strength / is 
shown in Fig. 5. Globally, the flow velocity Vf grows with 
increase of 7. Two regimes are clearly seen: Vf w ^(7) 
for (7) < 7 C and v f « (rjf (j) /m) 1 / 3 / 10 for (7) > j c « 
/ 1 / 2 /[30(r^m) 1 / 2 ]. In fact this dependence is very close 
to the one found in a model with fixed 7 |l6|. As a result, 
from (0) we obtain the dependence of flow velocity on 
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FIG. 5: (color online) Dependence of the absolute value of 
average velocity of particle flow Vf on average friction (7) 
for the parameters of Fig. 4 with / = 4, 8, 16, 32 (same 
symbols, from bottom to top). Dashed lines show the scaling 
dependence Vf ~ (/(t)) 1 ^ 3 at large (7) for different / values; 
the full line shows scaling Vf = (7) at small (7). Crosses 
show data for the same parameters as for squares (/ = 8) 
but with additional circular scatterer added in the center of 
unit cell to eliminate orbits with a straight flight through the 
whole system (see text). Insert shows the dependence of flow 
direction angle ip on polarization angle 0; data are given for 
/ = 8, uj = 1.5, r 2 = 50 and 4 < T < 11; full line shows 
average dependence <p = tt — 20. 

temperature: 

Vf/v w r d f/50T , (T<T C ); (4) 
Vf/v w (r rf //8T) 2 , (T>T C ); (5) 

where v = (2T/m) 1 ^ 2 is the thermal velocity and T c w 
Tdf is linked to j c obtained from Fig. 5. The transition 
between two regimes takes place when the energy given 
by ac-force to particle between two collisions becomes 
larger than thermal energy (T <T C ). In that case the ef- 
fect of driving is strong and Vf ~ fr c /m ~ r&f /(mT) 1 / 2 
leading to For T > T c thermal fluctuations are 

strong and the ratchet effect appears only in the sec- 
ond order of force / giving ((5J). The numerical factors in 
©,© are taken for the case = from Figs. 4,5. We 
note that Eqs. @-© are derived in the regime of rela- 
tively weak friction (7) <C uj and relaxation rate 1/r <C u. 
Another important point is that the dependence Q,© 
is robust with respect to variation of scatter geometry, 
e.g. introduction of additional disk scatterer in the cen- 
ter of unit cell eliminates all collisionless paths but gives 
no significant modifications (see Fig. 5). 

The dependence of flow directionality, determined 
through angle tp (vf = Vf (cos (p, sin y?)), on the polar- 
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ization of ac- force is shown in Fig. 5 (insert). In average, 
it is satisfactory described by the relation cp = tt — 26 
(similar dependence was seen in On a qualitative 

ground, we may say that at 6 = due to friction a particle 
becomes trapped between semi-disks of a unit cell that 
gives a directed transport to the left while for 6 = tt/2 
vertical oscillations push particle to the right in presence 
of friction. The linear dependence cp = it — 26 interpolates 
between these two limits. However, a more quantitative 
derivation is needed. 

It is interesting to apply the approach developed above 
to other types of thermostat. It is possible to realize 
the semi-disk Galton board with antidot superlattices 
for 2D electron gas in semiconductor heterostructures. 
With such structures the Galton board of disks has al- 
ready been implemented (see e.g. Il8l) and effects of 
microwave radiation has been studied . For disk anti- 
dots like those in 0, ^| the ratchet effect is absent due 
to symmetry of antidot. However, for semi-disk antidot 
lattice strong ratchet effect should appear. To find its 
properties we should take into account that in this case 
we have the Fermi-Dirac thermostat with the Fermi en- 
ergy Ep ^> T. Due to that in © the particle velocity v 
is equal to the Fermi velocity vp = (2Ep /m) 1//2 indepen- 
dent of T. This modification gives the average friction 
jf f° r the Fermi gas 

lF = Cfv F r d /T 2 w vj /r d , (6) 

where we kept the same numerical constant C ~ 1/50. 
In fact (0) follows from Dp ~ j 2 vprd (see (0) and 
^p ~ D E /T 2 . The second equality in © appears due 
to the fact that Ep ^> T c implying the regime Q with 
Vf ~ jF^d- Of course, only a small fraction T/Ep of elec- 
trons near Ep contributes to this ratchet flow. Hence, the 
current / per one semi-disk row is 

/ ~ er d n e v f T/E F ~ Cer 3 d ^T e f 2 /(Th) . (7) 

where we used that for the 2D electron Fermi gas Ep = 
7rn e h 2 /m. We note that in semiconductor antidot lat- 
tices like in 0, ^| the effective mass m is about 15 
times smaller compared to the electron mass. For a typ- 
ical parameters of semi-disk antidot lattice with electron 
density n e ~ 10 12 cm -2 , ~ 1/im, field strength per 
electron charge f/e ~ IV/ cm and T ~ 10 K we ob- 
tain Vf/vp ~ 10 -4 . At these parameters Ep ~ 150if, 
vf ~ 3 • 10 7 cm/ sec and the current / ~ 10~ 9 A is suf- 
ficiently large to be observed experimentally. The re- 
sult is based on the semiclassical estimate for the 
diffusion rate D E which assumes that the energy of mi- 
crowave photon is larger than the level spacing A inside 
one unit cell: Huj > A w 2ith 2 / (mr^) . In the opposite 
limit fajj <C A, ac-driving is in the quantum adiabatic 
regime when the excitation in energy is very weak. Thus 
for rd ~ ljim we have A w 5 • 10 _6 eV « 0.05i^ and 
the directed transport appears only for co/2tt > 1GHz. 



In experiments |(j the frequency was deeply in the adia- 
batic regime with uj/2tt ~ 100Hz and the directed trans- 
port was absent at zero mean force. We also note that in 
the quantum case the ratchet transport should disappear 
as soon as the amplitude of oscillations f /muo 2 induced 
by ac-force becomes smaller than the wavelenght h/mvp 
at the Fermi level. Thus the ratchet survives only for 
ou < y/ ' fvp/h. For field strenght of IV/ cm this gives an 
approximate bound at 30GHz. 

In summary, we showed that zero mean ac-force ap- 
plied to particles being in thermal equilibrium in asym- 
metric periodic potential creates a directed transport 
flow. Its direction is efficiently changed by polarization of 
the force. We also established the dependence of the flow 
velocity Vf on temperature and driving field strength for 
the Maxwell [Eqs. Q,©,®] and the Fermi-Dirac [Eqs. 
©,0] thermostats. 
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